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APPROXIMATE MODELS FOR CENTRAL SERVER SYSTEMS 



WITH TWO JOB TYPES 
by 

J. P. Lehoczky* 
D. P. Gaver** 



1 . INTRODUCTION 

1 . 1 General Summary 

The purpose of this paper is to provide a new approximation 
technique for describing the performance of a closed queueing 
network when there are two or more job types present, and the 
order of service (scheduling policy) of interest does not 
make system states Markovian in the number present at each 
server; "first-come, first-served" order is an example. The 
model also allows incorporation of non-exponential distributions 
to represent service time distributions. Our technique drastically 
reduces the size of the state space used to describe complex 
processes. Consequently, it facilitates calculation of 
state probabilities and system performance measures through use 
of numerical methods such as Gauss-Seidel iteration. It is also 
helpful in simulation studies. 

*J. P. Lehoczky acknowledges the support of the Air Force Office 
of Scientific Research at Carnegie-Mellon University, Grant 
No. AFOSR-74-2642B < 

* *D . P. Gaver acknowledges the support of the Command and 
Control Technical Center, Defense Communications Agency. 
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While we emphasize applications to models of multipro- 
gramming systems in this paper, the method described here 
may be employed for treating other job-shop-like problems. 



1 . 2 Background 

Multiprogramming computer systems have often been repre- 
sented by some variant of the classical Markovian cyclic queue- 
ing, or repairman, model described by Feller (1957). Early work 
by Gaver (1967) and Rice (1971) has been extended and improved 
upon by many others, e.g. Baskett et: al. (1975), Buzen (1973), 
Chandy (1975), Reiser and Kobayashi (1974) to name a few. 
Validation of these models, and applications to planning 
computer systems, have been reported by Giammo (1976), 
and Rose (1976). The models have, by and large, depicted jobs 
or programs as passing in a probabilistic (Markovian) manner 
from server to server in a closed network of queues forming 
before the various servers. The latter are identified as one or 
more Central Processing Units (CPU) , and one or more Peripheral 
Processors (PP) , such as disk or tape drives, or input-output 
devices. 

Papers by Jackson (1963) and Gordon and Newell (1967) 
early pointed out the simple "product form" structure of the 
long-run (or steady-state, or stationary) joint distribution of 
jobs present at each server in certain closed cyclic networks. 
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The latter solution holds if a fixed number of jobs remain in 

the system, and if job holding or processing times at Server i 

are identically and independently exponentially distributed, 

and transitions of jobs from server to server are Markov, 

each transition from Server i to Server j having the 

same probability P-^ j • F° r other conditions see Chandy (1975). 

The latter assumptions have frequently been made in modeling 

computer system behavior. Baskett, et^ a_l. have also obtained 

product-form solutions for situations involving different job 

types, so that the service time of Job Type k at Server i has 

rate parameter and the between-servers transition prob- 

(k) 

abilities are specified as p^^ . However, special queue 

disciplines must be assumed in order to achieve the results. 

The first-come, first-served (FCFS) discipline is not among these, 
except when identical service rates prevail. And yet FCFS is 
more representative of PP behavior than is, for example, a 
processor-sharing discipline. 

Computer system modeling based on direct application of 
the simplest Jackson-Gordon-Newell results is subject to question 
on several grounds. First, service or holding times need not 
be exponentially distributed (nor independent). Second, a 
mixture of job types may be expected in many real computing 
systems. For example, some jobs, e.g. those termed compute- 
bound , may require relatively extensive CPU activity, and spend 
relatively little time at the PP ' s ; others, the I/O or peripheral - 
bound, have characteristics that are the reverse. Also, in real 
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systems different jobs will make heavy demand for information on 

their own particular PP ' s , systematically ignoring others. This 

behavior is, in principle, not well-represented by the homogeneous 

character of the usual cyclic queueing assumptions. Additionally, 

non-exponential, e.g. hypo or hyperexponential, service times are 

of course not properly represented by simple exponentials. One 

remedy for this deficiency is to utilize phase-type distributions, 

see Cox and Miller (1965) and also Neuts (1977), One example of 

a phase type distribution is the gamma or Erlang, thought of as 

representing a sum of, say, k independently and exponentially 

distributed random variables. Specification of the phase of service 

in progress is all that is necessary to render the service 

th 

process Markov; when the k — phase terminates, service is complete. 
More generally, phase distributions arise as the hitting or first- 
passage time of a continuous-time finite Markov chain with one 
absorbing state; a service terminates when the absorbing state 
is hit. The class of phase distributions is large, containing 
distributions exhibiting features such as bimodality and nearly 
spiked (delta function) behavior. While Markov queueing network 
analysis to include phase type service is in principle not 
difficult, in practice, it may result in rather serious computa- 
tional problems, for the state space necessary to describe the 
system easily becomes enormous. This is especially true when 
it becomes necessary to include information concerning the order 
of jobs, classified as to types, in the state description. 

While such an "exact" analysis may be carried out for small 



4 



Markovian networks, it is advantageous to reduce the state 
space size in order to facilitate numerical calculations. For 
that purpose we therefore introduce and evaluate a weighted 
processor sharing FCFS model. 
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2 . A Numerical Approach Via an Approximating Process 



An appropriate Markov central server model of certain 
multiprogramming computer systems with realistic service protocols 
can be constructed by specifying an elaborate, large, state space. 

An example of such an approach is that of Gaver and Humfeld (1976) , 
and Humfeld (1977). In the latter investigations a computational 
method is devised for deriving the equations of probability balance 
for a system having two job types and FCFS discipline at one CPU 
and at each of k different PPs. There are J\ jobs (programs) 
of each type; i = 1,2. 

Such a model may, unfortunately, be far too large to be 
solved numerically for realistic systems. In the next section we 
introduce a stochastic process, Markov in the number of jobs of 
each type at each processor, that approximates the behavior of 
the more complex process required to describe the effect of FCFS 
disciplines. Apparently a product form analytical solution will 
not be available for our new approximating process. However, by 
eliminating reference to both number present and type order, the 
state space size is dramatically reduced. In terms of the 
approximating processes one may also more easily study many different 
queue disciplines as they effect performance in a multitype multi- 
programming environment. In effect, processor sharing, FCFS 
and various priority disciplines may be given a unified 
treatment. 
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3 . Weighted Processor Sharing FCFS Approximation 



For clarity of exposition, we assume that each of the 
service centers has a FCFS queue discipline and that there are 
two types of customers. The system contains a fixed number, 
say J\, of jobs of type i, i = 1,2. On device j, customers of 
type i have an exponential service time with parameter 
(j = 1 and 2 refer to PPl and PP2 , while j = 0 refers to the 
CPU). We define N^(t) to be the number of customers of type i 
at service center j, i = 1,2; j = 0,1,2. With a processor sharing 
queue discipline at device j , • { (N . . (t ) , N . _ ( t ) ) , t ^ 0 } 
forms a Markov chain where the probability a customer of type i 
completes service in [t,t + dt] is given to be 

y i N ji (t) dt / (N jl (t) + N j2 (t)) + o(dt) ' 1 = 1 / 2 • With a FCFS 
queue discipline, these probabilities are completely dependent 
upon the ordering in the queue. They are p-^dt and 0 if a 
type 1 job is in service, or 0 and if a type 2 job is in 

service. To preserve the Markov property and properly model the 
system one must keep track of the order of the jobs in the queue. 
This approach was adopted by Gaver and Humfeld (1976) for 
exponential service times. The resulting state space has 



J 1 + J 2 



J 1 + J 2 + k 



states if there are k distinct PP ’ s . This number must be 
greatly increased if service distributions of phase type are 
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allowed. This approach is thus seriously limited by the size 
of the state space. 

We would like to develop transition probabilities which 
approximate FCFS behavior but which are based only on N_.^(t), 
i = 1,2; j = 0,1,2, not the ordering in the queue. Let us focus 
on a single service center with 1SL (t) jobs of type i present 
at time i and exponential (y^) service rates (i = 1,2). In 
the spirit of processor sharing, we approximate FCFS as follows: 
let A^(t) = N^(t + dt) ” lSL(t), i = 1,2. We assume 



P(A 1 (t)=-l, A 2 (t)=0 |N(t) ) 



N 1 (t) /v 1 
N" 1 (t)/y 1 + N 2 (t)/y 2 



dt + o(dt) 



P ( A x (t ) =0 , 



A 2 (t)=-l|N(t) ) = y 2 



N 2 (t)/y 2 

N 1 (t)/y 1 + N 2 (t)/y 2 



dt + o(dt) 



P(A 1 (t)=0, A 2 (t)=0 | n (t) ) 



N x (t) + N 2 (t) 
N 1 (t)/y 1 + N 2 (t)/y 2 



dt + o(dt) . 



(3.1) 



The motivation for the suggested approximation is as 
follows. Over a short period of time [t, t + dt] , only the job 
currently in service has positive probability of completing 
service. This is in sharp contrast to (3.1) where both types have 
positive probability of a service completion. Nevertheless, if 
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we focus on a longer period of time, say the amount of time 

needed to service all customers currently in the queue, then on 

the average N^(t)/y^ time units will be spent servicing type i 

jobs, i = 1,2. Hence the server will spend approximately a 

fraction (N. (t) /y . ) / (N, (t) /y. + N_(t)/y 0 ) of this total time 
1 11 1 z z 

servicing type i customers. Approximating the system using (3.1) 
amounts to applying this fraction, which changes dynamically, to 
every interval of length dt , rather than just to the longer 
interval in which every customer is served. The approximation 
is intended to be applied only for steady state calculations, 
since it does not model the actual dynamics of a FCFS queue 
discipline faithfully. We note in passing that this is only a 
first order approximation. It may be improved by applying the 
fraction E(X^/(X^ + X^) ) to each dt time interval where X^ 
and X 2 are independent gamma random variables with parameters 
(N^(t),y^) and (N^Ct)^^). This tends to improve the approxi- 
mations, but at the cost of a much more complicated analysis. It 
may, however, be quite important when phase type service distributions 
are assumed. 

Equations (3.1) are a special case of more general tran- 
sition probabilities for single server queues with multiple 
customer types. Suppose there are T types and (t) is the 
number of type i customers (1 <_ i <_ T) enqueued. Let A^(t) 

= ISL(t + dt) - ISL(t). Equations (3.1) can be generalized as 
follows 
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P(A i (t)=-l, Aj(t)=0, j / i | N ( t) ) 



dt + o(dt) 



N i (t) 

y T . a. .N . (t) 
£ D=1 ID D 



(3.2) 



P(A i (t) = 0, 1 < i < T|N(t) ) 



= 1 — dt + o (dt) , 



y T 

1 j=i 



a. .N . (t) 
i] ] 



Equations (3.2) provide a parametrized family of single 
server queues. By adjusting the parameters { a.^ , 1 £ j <_ T} 
one can create many different service disciplines. For example, 
if a^j = the resulting discipline is approximately FCFS. 

The case a^ = 1 corresponds to processor sharing. Typically 
one would choose a. . = +1 and a. . = 1/a. . . By letting a. . 

li ij ji J 3 ij 

approach 0 or 00 a priority (of i over j or vice versa) 
will be created; as a^ -* 0, type i jobs are given priority 
over type j . The coefficients a^ ^ provide relative weights 
for each type of job at a particular server. See previous work 
by Lehoczky and Gaver (1977) utilizing a similar approach; that 
paper also presents numerical validations. 

The parameters a^ have intuitive meaning as explained 
before, and the parametric framework allows for the possibility 
of optimizing some system measure of effectiveness (such as 
idleness or queue lengths) by appropriate choice of the a ^j* 
Once an optimal set of a^^ has been found, it can be converted 
into a scheduling algorithm which allocates a particular weight 
to each type of customer. 
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This methodology (equations (3.1) or (3.2)) can easily 
accommodate phase distributions, but at the price of a possibly 
great increase in the size of the state space. Suppose that at 
a particular server the service distribution for type i customers 
i of phase type with n^+1 states, initial distribution a^, 
service or holding time vector y^ and transition probabilities 

, 1 < i < T. We identify each state of each phase distribution 

rT 

as a type of job: there are S = z,^_^ ru types. Let ^ (t) 

be the number of type i jobs in phase j at time t, 

n . 

1 _< j _< n^, 1 £ i £ T. Here \ , (t) = N^(t), the number 

of jobs of type i. 

Coefficients {a.. 1 _< i , j _< S) can now be defined by 

establishing the relationship between types i and j. Equations 
(3.2) can be easily applied with T replaced by S. While the 
number of states is potentially enormous, many interesting cases 
can be handled by solving the steady state balance equations by 
Gauss-Seidel iteration methods. 

We have described the situation of a single service center. 
The central server model consists of many servers (CPU and PP's). 
The servers are assumed to be linked in a Markov fashion. When 
a type i job completes CPU service, it moves to PP . with prob- 
ability a ^j / a ^j £ 0, a ^j = + 1 • Upon completing service 

at PP j , it moves back to the CPU. The number of jobs of each type 
is typically held fixed, but this can be modified to allow them 
to fluctuate, with the total fixed, in a straightforward fashion. 
One can also allow the total number of jobs in the system to 
vary with time although such models are of lesser importance. 
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4. An Illustration 



To illustrate the quality of the approximation given by 
(3.1) we give an exact analysis of a very simple situation in 
which closed form expressions can be obtained. A detailed numerical 
analysis of a more realistically sized system will be given sub- 
sequently . 

We study the simplest possible multitype central server 
model with a CPU, two PP ' s and two job types, each with one job. 

To further simplify the calculations, type i jobs always move 
to PPi, i = 1,2 and then back to the CPU. Assume service rates 
A^ , i = 1,2 at the CPU and 1 at each PP . 

Assuming the CPU is operating in a FCFS fashion, the state 
space has 5 states defined by the number of each type at the CPU 
and the type of job currently in service if both types are present. 
Here x^ will be used to denote the long-run probability that 
the system inhabits the state having label i. 



STATE 


LABEL 




BALANCE 


EQUATIONS 


(0,0) 


1 




2X 1 


= 


A r 


X 2 


+ X 2 x 3 


(1,0) 


2 


(x i 


+ 1) x 2 


= 


X 1 


+ 


X 2 x 5 


(0,1) 


3 


u 2 


+ l)x 3 


= 


X 1 


+ 


X l x 4 


(1,1,1) 


4 




X l x 4 


= 


X 2 






(1,1,2) 


5 




X 2 x 5 


= 


X 3 












1 


= 


X 1 


+ 


x 2 + x 3 



(4.1) 
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Equations (4.1) can be routinely solved to give 



x = (a 1 a 2 + x 1 + a 2 )/d , 

x 2 = (2 + A 2 ) /D , 
x 3 = (2 + A 1 )/D , 

x 4 = (2 + A 2 )/A 1 D , (4 ‘ 

x 5 = (2 + A 1 )/A 2 D 2 , 

D = [(A 1 + 1 ) 2 ( A 2 + l) 2 - 1]/A 1 A 2 

The model has only 4 states if the approximation outlined 
in (3.1) is used. States 4 and 5 in (4.1) merge into a single 
state; this represents the advertised economy of the proposed 
approximation. Here y^ denotes the long-run probability of 
inhabiting label state i. 

STATE LABEL BALANCE EQUATIONS 



(0,0) 


1 


( A-^ + 


2y l = X l y 2 


+ A 2 y 3 


, -1 


(1,0) 


2 


+ 

I 1 

>1 

II 

CM 

>1 

i — 1 


(1/A 1 + 


1/A 2 ) y 


(0,1) 


3 


(a 2 + 


i)y 3 = y 3 + 


Y 4 <1/*! 


+ i/a 2 )“ 


(1,1) 


4 


2(1/A 1 + 1/A 2 ) 


" ly 4 “ y 2 + 

1 = y-L + 


y 3 

y 2 + y 3 


+ y 4 



(4.3) 
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The solution is easily found to be given by 



y-j_ — ^j_^2^ 2X 1 X 2"^" X 1 X 2 ^ 

y 2 = 2X 1 X 2 (X 2 + 1)/E 

y 3 = 2X 1 X 2 (X 1 + 1)/E (4.4) 

y 4 = (X 1 + X 2 } (2 + X i + X 2 )/E 
E = ( 1 ^X 2 + A ^ + A 2 + 2 ) ( 2 X^X 2 + A^ + X 2 ) 

One can compare the results by aggregating and x,- , 

so the state spaces are the same. First x^ = y^, the state which 
corresponds to CPU idleness. For the other three states the 
approximation is not exact. We give the percent error of y^ 
computed by |x^ - YjJ/Y^ where x^ is the aggregation of x^ 
and x,- . 



|x 2 - y 2 |A 2 = a 2 | x 2 - x 1 |/(2 + x 2 ) (2x 1 x 2 + x x + x 2 )) 




|x 4 -y 4 |/x 4 = X 1 X 2 (X 2 - X 1 ) 2 /(2X 1 X 2 + X x + X 2 ) (X 2 + 2X 1 + X 2 + 2X 2 )) 

(4.5) 

One may make a detailed analysis of (4.5). The overall 
result is that the percentage errors are very small unless an 
extreme case is chosen. To illustrate we focus on |x 2 - y 2 l/x 2 
and let X 2 = kA^. In this case the percent error is given by 
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(4.6) 



|x 2 - y 2 | k(k-l)A 1 

x^ = (2 + k X' ) (2k A ^ + k + 1) 



For fixed k, this expression can be maximized as a function of 
A^ by selecting A-^ = /k + 1/k . Substituting this value into 
(4.6) gives 



|x 2 - y 2 



-max = 



k - 1 



(2 + /k+1) 



(4.7) 



The following table shows the maximum percent error as a function 
of k. 



_k 

1 

2 

3 

4 

5 

10 

50 

100 



(k-1 ) / (2 + /k+T) 2 
0 % 

7.1% 

12.5% 

16.7% 

25.2% 

31.8% 

59.8% 

68.9% 

100 % 
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One can see that except in very extreme cases, the percent 
error will be very small indeed. A similar analysis can be 
carried out on the other terms with similar results. 
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5 . A Second Approximation 

Another approximation method can be used which offers 
several advantages, but some disadvantages, over the weighted 
processor sharing example. We illustrate this for a single service 
center. Rather than keep track only of the number of each job 
type at the service center as in the first approximation, we expand 
the state description to include the number of each type awaiting 
service and the identity of the customer in service. Arrivals 
to the queue and departures from the queue are easily handled 
since the identity of the customer in service and the number at 
the center is known. An approximation is required to identify 
the type of the next customer to enter service when a customer 
leaves, because the order within the queue is not known. We 
adopt a random service order policy approximation. Specifically, 
if W^(t) represents the number of type i jobs waiting for 
service, i =1,2 and I(t) = i if a type i job is in service, 
then we assume, letting A^ = W. (t + dt) - W^ (t) 



P(A 1 =-1 



V 0' 1 (t+h) =1 |w(t) , X(t)-l) ^ u iM+ „ 2(i) 



Wj^t) 



dt + o(dt) 



P(A 1 =0, A 2 =-l, I (t+h) =2 | W (t) , I (t) =1) = p 




dt + o(dt) 



1 w x (t)+w 2 (t) 



W, (t) 



P(A 1 =-1 



A 2 0, I (t+h) l|W(t), I(t) 2) p 2 ( t ) +w 2 (t) 



dt + o(dt) 



P(A 1 =0, A 2 — 1, I (t+h) =2 | W (t ) , I (t ) =2 ) = p 



w 2 (t) 



dt + o(dt) . 



2 W x (t)+W 2 (t) 



(5.1) 
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This can be easily extended to many service centers and provides 
a second approximation to FCFS service disciplines. The numerical 
accuracy will be assessed in the next section. 

This approximation requires a slightly larger state space 
than the weighted processor sharing (3.1) approximation for 
exponential service distributions. Nevertheless if phase type 
distributions are introduced the state space for (5.1) is much 
smaller than for (3.1). Phase type distributions can be very 
easily accommodated. Another advantage is that (5.1) provides 
better accuracy than (3.1), although both give excellent results. 
On the negative side (5.1) cannot be conveniently changed to 
allow for other types pf queue disciplines. 
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6. Numerical Illustrations 



Numerical studies were carried out using (3.1) and (5.1) 
and were compared with the results of Gaver and Humfeld (1976) . 

Below we present the results for three cases. In all cases there 
are two PP ' s and two job types. The CPU service rates are A ^ 

while (jj^, are those for PP1 and (y 2 i ' M 2 2 ^ for pp2 • 

and a.^ are the probabilities a type 1 or type 2 job moves to 
PP1 upon leaving the CPU. is the fixed number of type i jobs, 

i = 1,2. 

Method 1 refers to the results of Gaver and Humfeld (1976) 
obtained by a Gauss-Seidel iteration of the balance equations for 
the exact system. Method 2 refers to the steady state results 
obtained using a Gauss-Seidel iteration of the balance equations 
derived from (3.1) , the weighted processor sharing approximation. 
Method 3 refers to the approximation based on a random source 
order (5.1). The results were obtained using simulation. The 
results are based on 120,000 system transitions, and estimated 
standard errors are provided with each quantity. 

The numerical values given in the previous tables illustrate 
the exceptional accuracy of both approximation method across a broad 
number of cases. The largest absolute error is .02, thus the 
approximations are extremely useful. The tables illustrate that 
method 3 (based on (5.1)) is more accurate than method 2 (based on 
(3.1)). Method 2 allocates slightly too much time to slow jobs, 
slightly too much to fast jobs. This results in the slight but 
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systematic differences between methods 1 and 2. In view of 
the minor differences, the complicated correction mentioned in 
section 3 is certainly not called for. 

The tables illustrate that method 3 should be used when- 
ever possible. This will be especially true if phase type dis- 
tributions are required. The accuracy is better than method 2 
and the state space requirements may be substantially reduced. 
Method 2 can, however, serve as a starting place for calculating 
approximate closed form expressions for occupancies and idleness 
probabilities, say by using a diffusion approximation analysis. 
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Occupancy 




Idleness 




( V J 2> 


Method 


CPU 


PPl 


PP2 


CPU 


PP1 


PP2 


(1,D 


1 


.830 


.582 


.588 


.406 


.519 


.501 




2 


.827 


.566 


.607 


.403 


.524 


.483 




3 


. 830 ( .003) 


.587 (.004) 


. 587 ( .006) 


.409 ( .001) 


.518 (.002) 


. 505 ( .005) 


(1,2) 


1 


1.172 


.613 


1.215 


.336 


.537 


.284 




2 


1.163 


.590 


1.247 


.334 


.544 


.266 




3 


1 . 172 ( . 013) 


. 616 ( .003) 


1 . 219 ( . 013) 


.337 ( .007) 


.535 (.004) 


. 283 ( .001) 


(2,1) 


i 


1.342 


1.093 


.565 


.282 


.341 


.556 




2 


1.343 


1.074 


.582 


.277 


.345 


.538 




3 


1. 346 ( .012) 


1.085 (.019) 


.576 ( .007) 


.281 (.006) 


.342 (.005) 


.533 ( .005) 


(1,3) 


1 


1.444 


.613 


1.943 


.305 


.558 


.165 




2 


1.425 


.587 


1.988 


.306 


.565 


.151 




3 


1.444 ( .019) 


.627 ( .006) 


1.940 (.018) 


.309 (.006) 


.556 ( .004) 


.170 ( .003) 


(2,2) 


1 


1.779 


1.106 


1.115 


.224 


.367 


.355 




2 


1.777 


1.077 


1.146 


.220 


.372 


.334 




3 


1 . 794 ( . 024) 


1.104 (.010) 


1.112 (.015) 


.221 ( . 003) 


. 369 ( .005) 


. 359 ( . 005) 


(3,1) 


1 


1.839 


1.626 


.535 


.218 


.245 


.594 




2 


1.844 


1.608 


.547 


.214 


.248 


.579 




3 


1. 826 ( .009) 


1.637 ( .009) 


. 547 ( . 009) 


.221 ( .002) 


. 242 ( . 002 ) 


.591 ( .003) 


(3,3) 


1 


2.847 


1.572 


1.581 


.138 


.296 


.286 




2 


2.846 


1.536 


1.618 


.134 


.300 


.268 




3 


2 . 836 ( .076) 


1.600 (.036) 


1.583 ( .042) 


.140 (.009) 


. 297 ( . 006 ) 


.292 ( .008) 


(2,5) 


1 


2.672 


.998 


3.330 


.182 


.441 


.096 




2 


2.619 


.959 


3.422 


.182 


.447 


.085 




3 


2 . 684 ( .048) 


1 .021 ( .022 ) 


3.315 ( . 047) 


.180 ( .006) 


. 438 ( . 003) 


.100 (.006) 


(5,2) 


1 


3.459 


2.690 


.851 


.115 


.169 


.491 




2 


3.475 


2.664 


.861 


• 111 


.172 


.477 




3 


3.456( .078) 


2 . 709 ( . 077) 


. 854 ( . 025) 


. 118 ( .005) 


.172 (.008) 


. 493 ( . 014) 
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Case II 


X 1 


= 0,526 




\ 2 = 0.339 


a ll 


= 1.0 






^11 


= 1.0 




y 12 = 1.0 


P 

to 


= 0.0 






y 21 


= 1.0 




y 22 = 1,0 










Method 


Occupancy 
CPU PP1 


PP2 


Idleness 
CPU PP1 


PP2 


(1,1) 


1 


1.610 


.200 


.190 


.059 


.800 


.810 




2 


1.608 


.205 


.187 


.059 


.795 


.813 




3 


1.609 ( .005) 


.20K.002) 


.192 (.003) 


.060 ( . 002) 


.800 ( .002) 


. 809 ( .0( 


(1,2) 


1 


2.575 


.135 


.290 


.014 


.865 


.753 




2 


2.571 


.136 


.292 


.014 


.864 


.754 




3 


2 . 574 ( . 001) 


.135(.003) 


.293 ( . 002) 


.014 ( .001) 


.866 (.003) 


. 751 ( . 00 


(2,1) 


1 


2.500 


.353 


.148 


.018 


.713 


.852 




2 


2.503 


.351 


.146 


.018 


.709 


.854 




3 


2.503(.011) 


. 351 ( .010) 


.149 ( .002) 


.018 ( .001) 


. 715 ( .006) 


.852 (.00 



Case III 





*1 


= 0.6 




X 2 = 0.8 


tt ll 


= .55 






y ll 


= 1.2 




Pl 2 = 1.0 


a 21 


= .30 






y 21 


= 0.9 




y 22 = 0.5 














Occupancy 




Idleness 






Method 


CPU 


ppi 


PP2 


CPU 


PPI 


PP2 


(4,4) 


1 


6.107 


.345 


1.548 


.oil 


.742 


.372 




2 


6.152 


.347 


1.501 


.009 


.742 


.376 




3 


6.113( .062) 


.342 (.015) 


1.561 (.060) 


.012 (.003) 


. 744 ( .008) 


.374 (.00 
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